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A DYNAMIC EQUATION FOR STOCHASTIC 
MAGNETIC FIELD LINES IN THE GALAXY 

Frank C. Jones 
Theoretical Studies Branch 

ABSTRACT 

It has been pointed out by previous authors that the turbulent motion of the ion- 
ized interstellar gas wiU cause a given line of force of the galactic magnetic 
field to wander in a random manner even though the average field may be quite 
regular. We derive here an equation that describes the development of the 
probability density of a given field line as a function of distance along the aver- 
age field. We demonstrate that a Fokker- Planck type of equation is quite un- 
suitable in this case and that the answer must be sought in the theory of station- 
ary Gaussian processes. We also derive the distribution function for the "end 
points" of a field line that starts on the galactic central plane. End points in 
this context means the point where the field line first wanders so near the sur- 
face of the galactic disk that it can no longer confine cosmic- rays. 
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A DYNAMIC EQUATION FOR STOCHASTIC 


MAGNETIC FIELD LINES IN THE GALAXY 
I. INTRODUCTION 

The earth is situated near (~ 20 pc) the central plane of the galactic disk. The 
average galactic field (~ 2-3 /> gauss) appears to lie mainly in the disk and run 
approximately parallel to the spiral arms (Davis and Greenstein 1951; Hiltner 
1956; van de Hulst 1967; Davis and Berge 1968). There are also theoretical 
arguments that it should lie essentially in the azimuthal direction (Parker 1969, 
Roberts and Yuan 1970). However, a given flux tube such as the one that passes 
through the earth and its immediate vicinity will not, in general, follow the pat- 
tern of the average field but will wander in a random manner about the galactic 
disk. This random wandering of a given Reid line or group of lines has been 
discussed at length by Jokipii and Parker (1969 a, b) and has been related by 
them to the problem of cosmic- ray propagation in tue galaxy. Briefly the re- 
lationship is as foUows: the random wondering of the field line to the surface of 
the galactic disk provides a means whereby a cosmic- ray particle, which at the 
position of the earth is buried deeply (~10 7 cyclotron radii) in the disk, may es- 
cape the disk in a time of the order of 10 6 years. At the point where the Reid line 
reaches the "surface" the cosmic- ray particles can escape presumably via the 
bubble blowing instability (Parker 1965). Since the galaxy does not have sharp 
edges the notion of the surface of the disk is necessarily a bit vague but it is 
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usually assumed that when a field line reaches a point sufficiently far from the 
central plane ( ~ 100 pc) the instability will occur and the cosmic- ray particles 
will escape. 

The positions of the ends, for cosmic- ray purposes, of the flux tube passing 
through the earth are important parameters in any discussion of the cosmic- 
ray anisotropy, whether from a smoothed out source point of view (Kulsrud and 
Psarce 1969) or from a statistical source point of view (Jones 1970 a, b). In 
either case the average cosmic- ray streaming observed at the earth will depend 
on the position of the earth relative to the end points of its associated flux tube. 

Unfortunately a determination of the actual positions of these end points is im- 
possible, at least by present techniques, and we must turn our attention to their 
probability distribution function. Jokipii and Parker (1969 b) discuss this ques- 
tion and give a rough estimate of where a given field line might end. Their esti- 
mate was very crude, however, and was of no value in determining where the 
maximum likelihood position would be. 

hi this paper we will consider this end point problem in more detail and obtain 
a probability distribution function for the end points of the flux tube passing 
through the earth. 

ha their analysis Jokipii and Parker (1969 b) determine that a field line whose 
mean position lies in the central plane of the galactic disk has a probability of 
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being found between z and z + dz above (or below) the central plane given 
by 

f (z) dz = exp (- ! z i z 0 ) dz z 0 . 

The general problem we wish to consider is: given that the field line in question 
is at z - 0 at a point y - 0, where y measures the distance along the mean field, 
what is the probability distribution in z as a function of y. Our specific question 
is the following: given that the field line is at z - 0, y = 0, the position of the 
earth, what is the probability that it continues a distance y down the mean field 
direction without reaching a distance ! z ! - z } from the central plane and thus 
terminating its role as a cosmic-ray container. 

The equation that we shall derive will be a dynamic equation in the sense that 
it will describe the development of the probability function along the direction y 
starting with some given form at y = 0 in a manner completely analogous to the 
development in time of dynamic stochastic systems. We shall not consider the 
actual time development of the stochastic field in this paper. Typical rms 
velocities of the gas clouds that carry the magnetic fields are ~ 5 Km/sec which 
correspond to motions of about 15 pc in a typical cosmic ray lifetime of 10 6 
years. We cannot claim that such a motion is truly negligible, however, we would 
not expect a pronounced effect from ignoring it at the present time. 

We shall first, following Jokipii and Parker (1969 a), consider the dynamic equa- 
tion to be of the Fokker- Planck type. The requirements that the asymptotic 
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probability distribution be of the form exp (-lzi/z 0 ) requires that the diffusion 
coefficient be of the form exp (lz!/z 0 ) (indeed any asymptotic form d(z) with the 
requirement of no probability current requites a diffusive coefficient of the form 
d~ l (z)). it will turn out that because of this and the fundamental assumptions in* 
herent in a Fokker- Planck type equation the results will be quite unsatisfactory. 
They will indicate a far too rapid transport of the field lines to the boundary. At- 
tempts to patch up this problem will be seen to be also unsatisfactory. 

It will turn out that we can construct a model of the wandering of field lines based 
on the theory of Gaussian processes from which we can extract a dynamic equa- 
tion that is free of the objectionable characteristics of the Fokker- Planck approach. 

II. THE FOKKER- PLANCK APPROACH 
The probability that a given field line can be found between z and z + dz above 
the galactic central plane at a point y along the direction of the mean field will 
be given by f(z; y) dz where f is the probability density function! We shall assume 
initially that f satisfies a Fokker- Planck equation of the form. 

1 

We shall adopt the notation of Feller (1966) in distinguishing between o probability density f(x) 
of a random variable, 

P(0 £ X - b) = I f(x)dx, 

•'a 

ond it* distribution function F(x), 

P(X < x) • F(x) = f f (y) dy . 

The various variables of a multivariate density or distribution function will usually be separated 
by commas though we will sometimes use a semicolon to indicate that the variables to its right 
are to be considered as parameters in this particular instance. In a conditional probability density 
the variables to the right of a vertical slash are the one that are fixed, i.e., 

P(x <X<x + dx iY = y) = f (x I y ) . 
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_L (z ‘> > + ~~ (d (Z)f(z.y))--^ — -J (d 2 (z)f (z . y)) 
«Jv *z 1 2 ‘z^ 


0 


( 1 ) 


where d t (z) & z /by and d 2 (z) «<&**> /• y. We require that 


f (* . y ) - exp (- i z i z 0 ) /2z o as y 


we therefore have 


[d 1 (z)exp(-|z| z 0 J- ~ - [d 2 (z)exp(-|z| /z 0 )] - 0. 


2 dz 2 


(2) 


The first integral of equation (2) is the probability current or 


d.(z)exp(-|z| z 0 ) ~ [d (z)exp(-iz|/z 0 )J = j = const. (3) 

* Z oz * 


j could be an overall constant which would destroy the symmetry of the problem 
with respect to plus and minus z and would add a decidedly unphysical aspect to 
the solution. A i which changed sign at z = 0 in order to preserve the symmetry 
would involve a source or sink of probability at z = 0 and would be equally un- 


physical. We therefore choose j = 0 everywhere. 


! , 


From (3) with j = 0 we obtain finally. 




( 4 ) 


where 


d 2 = d(d 2 >/dz and e(z) - ±\ for z < 0 

respectively. If we divide the galactic magnetic field into an average field B A 
and a small fluctuating part B , we have ( Jokipii and Parker 1969 a) 

d j (z) = <bz>/by - <B,>/B 0 
d 2 (z) = <Sz 2 >/5y £ 2L <B? f >/Bj 

where L is the correlation length of the random field B , . Since we must have 

<V • B,> = V • <B.> = 0. 

this gives 


d<B lt > B<B ly > 
3z 


(we are considering two dimensions only but similar considerations would apply 
to the three dimensional case). If we wish our problem to be stationary, i.e. all 
statistical parameters are independent of y we must require 


3<B > 3^B. > 

LL. 12L. = 0 and < B . > = con s t . 

3z 3y 
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Once again we require * B,> * 0 to preserve symmetry and we are left with 


d 2 (z) 


— d 2 (z) - d 2 (z) - v exp( |z | z 0 ) 


where 


v - (< i*>z 2 > /by ] f a0 


We have arrived at the equation 


3f(z.y KZJL f exp( | z l / ? 0 ) f ( z . y ) ] -- 0 
Sz 


2 : 4 ,2 


( 6 ) 


with the initial condition f(z, 0) - <? (z) as our fundamental equation. If we make 


the substitution 


f(y, z) = exp(-| z j /z 0 ) u(z.y) 


we have the equation 


exp (— 1 z I /z 0 ) 0 

0y 2 dz 2 


( 7 ) 


with the initial condition 


u(z) 0 - 6(Z) 


and asymptotic condition 


Equation (7) may be solved by the standard techniques of constructing a Green 
function after Laplace transforming the y variable. The solution is 


u( I * I .y ) 


1 

2z, 


L 

n ^ 0 




exp(- \z\Zz 


»>J 


Jo 


exp(- 


: 2 y 


8z‘) 


<«> 


where J Q (x) is the Bessel function of order zero that is regular at x « 0 and 0 n 
is the nth root of the equation (x) « 0 with p 0 - 0. It is easy to see that the 
effective length for relaxation of the initial probability distribution is y 0 - 
8 z * 

We now wish to use the dynamic equation (6) to compute the probability F (y; z j ) 
that a field line has reached the point y without ever having made an excursion 
in | z | as great as z . If one thinks a moment about standard particle diffusion 
with totally absorbing walls one will realize that after a time t all those particles 
that are left are just those that have survived a time t without ever having 
reached the absorbing wall. We therefore wish a solution of equation (6) with the 
boundary condition 


f , (z , y ) = 0 for z = ±Zj 


FjCyiZj) = 



f i(z,y) dz. 


and then 



The distribution function for the endpoints is 

F^fy; z ,) - 1 - F,(y ; z t ) 
and since the endpoint density is 

f 2 (y • z i ) " J F 2 /^y 


we have using (6) 


f 2 <y= z i> = -a7 F .(y ; u> = - j — — d * 




H 


d 2 


dz' 


[exp( I z | /z 0 ) f ,(z.y )] dz 


- Z , 


- V I du 

" 2 V dz 


c*u 

dz 


= v 


du 


djz! 




Once again using standard techniques we find that f (y; z t ) is given by 


do 


f 2 (y ; z i> = 


V 


L 


exp [-vB 2 y/8*J] 


2rz l 77 1 CdD/dnl^' 


where the p are the zeroes of the function 


D(^; 2,) s J,(*7) Y 0 (i7x t ) - Y t (77) JqOjXj) 


where Y is the mth order Bessel function of the second kind and x . a 



One can now compute the distribution function of the end points i.e., the proba- 
bility that the position of endpoint Y is less than or equal to y: 


P(Y < y ) = F 2 (y; z t ) 


However, since this approach will turn out to be unsatisfactory we will not waste 
time examining F 2 (y; z x ) in detail but will take a simpler approach. 

It is easy to see that the various moments < (y/y 0 ) m > are given by 


f 


(y/y 0 ) m> s (y/y 0 ) m f 2 ^ ; 2 i) d Y 


(ii) 


1 m ! T — ; 

77 trl (V ) m+1 c 


[dD/d rt] Vs 0 


and they are, of course, functions of z x /z 0 . 


In Table I we have listed various values of z { /z Q and the corresponding values 
of < y/y 0 > or the mean length of a field line. 

One notices immediately that the mean lengths of a field line are considerably 
less than the characteristic length y 0 even for an absorbing wall as much as 
four characteristic heights z Q away. To evaluate the length y 0 = 8zJ/v we note 
that 


*;s2L<Bj t >/B2| t ^ But <B| i >/Bq 


*0 




z*/V 


10 


80 


y 0 ft 2L 

This says that for z l /z 0 - 2.2 (the case considered by Jokipii and Parker 
1969 b) an average field line will reach one of the surfaces z = ±z x in a distance 
< y > = .65L. 

Jokipii and Parker's estimate was about ten times this or 6L. Furthermore 
this says that an average field line has a slope 2.2 z Q /.65L $ 3.4 z Q /L whereas 
we initially assumed a mean slope of ft z^/L so we have an internal contradiction 
as well as an external one. 

The difficulty arises because we have a situation that violates the basic assump- 
tion underlying all diffusion type equations. The assumption is that there exist 
at least two scales in the problem: a microscopic scale of the forces driving the 
system (the random field B is the force in our problem), and a macroscopic 
scale on which the forces and the response of the system are completely un- 
correlated. The diffusion type equation describes the behavior of the system on 
the macr oscopic scale only. We can see that this is not the case here; the 
microscopic and macroscopic scales are both characterized by the correlation 
length L. It is this violation of the two scale requirement that causes the Fokker- 
Planck approach to yield self contradictory results. 


One could attempt to patch up the Fokker- Planck equation by replacing the 
constant diffusion coefficient i • = 2L < B} 2 > / B 2 by a generalization that is y 
dependent, 


t '(y) = (2'B 0 2 ) C ll (y')dy‘ 

•*0 

where C zz (y) is the zz component of the random field correlation tensor, i.e. 

C ij (y ) = <B i(x) B. (x + y) > 

Using v (y) in equation (1) is equivalent to replacing the variable y with the vari- 
able v where ??(v) = (^ (<*))” 1 1 ^(y ') dy' in the original solution of equation (1). 

•'o 

This is only an ad hoc procedure and cannot be considered completely satisfac- 
tory. It does, however, modify the rapid initial dispersion of probability, the 
dispersion is proportional to y rather than y 1/2 for small y. The ove all effect 
is small though and only increases the mean field line length by an addition of 
approximately L/2. 

We should, therefore, turn to an entirely different approach to the problem, one 
in which there is no need to postulate two length scales, a microscopic scale and 
a macroscopic scale. Such a procedure can be found in the theory of stationary, 
Gaussian processes which we shall discuss in the next section. 
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III. THE GAUSSIAN PROCESS APPROACH 


The random wandering of the galactic magnetic field is produced by the turbulent 
motion of the gas to which the field is frozen. We shall, therefore, turn our at- 
tention to a particular model of this underlying interstellar gas. This model 
will not necessarily be a correct picture of the turbulent galactic gas, however, 
it is at least a physically realizable model and from it we can calculate, without 
any mathematical approximations, the statistical properties of the magnetic Held. 
This means that while the results we obtain for the field line wandering may not 
correspond exactly with reality, they will at least be physically realizable and 
will contain no unphysical properties introduced via some mathematical approxi- 
mation scheme. As we have seen in the last section not all methods of generating 
dynamic equations for probability densities can make this claim. The hope is, 
of course, that our results will constitute a reasonable approximation to the situ- 
ation that exists in our galaxy. 

Consider the following two dimensional situation. A quiescent gas has a uniform 
magnetic field Bey embedded in it. The field lines may be labeled by the value 
of their z coordinate at the y = 0 plane. Now imagine that the gas is displaced 
in the z direction by an amount s which varies from one point to another depend- 
ing on the value of the y coordinate thus s ■ s(y). We now consider an ensemble 
of such situations and let s(y) be a random function on this ensemble. We shall 
assume s(y) to have the following properties: it is expressible as a Fourier 
series 
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s(y) 



S t cos (kj y + (p x ). 


and the phase angles 4> i are random variables, uniformly distributed between 
0 and 2n, From these properties it can be shown (Rice, 1944, 1945) that s(y) is 
a stationary, Gaussian processes. This means that the multivariate probability 
density for the N quantities x Jt x . . . x N where the x . are values of the 
function s or any of its derivatives with respect to y taken at N points y , 
y . . . y (where some or all of the y. may be identical) is given by the N 
dimensional Gaussian density 


f (x) dx - (2 7r)" N/2 |mI ” ,/2 exp 

'V •'V 





( 12 ) 


where x and x are the row and column vectors respectively with components 
x , x ... x and M is the matrix M = < x x >. The process is stationary 

12N i j i j 

because the quantities < x . x > are functions of the relative separation only. 
From this it is evident that the probability density function is a Gaussian rather 
than an exponential as inferred by Jokipii and Parker (1969 b). Of course, what 
is observed is the velocity distribution of the interstellar gas clouds which is 
approximately an exponential function of cloud velocity (Spitzer 1968). The 
position and velocity distributions will have the same form only if the dynamic 
(in the time development sense) system is a linear one and it is tempting to 
argue that non-linearities could cause the velocity distribution to be exponential 
even though the position distribution is Gaussian. 
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This argument is dange^uus, however, for if the system is appreciably non-linear 
the Fourier modes will not be normal modes in the dynamical sense and even if 
they are started with random phases mode coupling can produce correlations in 
the phases as time advances. So it appears that we must simply acknowledge 
that this model departs somewhat from observation and hope that the departure 
ia not too serious. 

From equation 12 the distribution of probability for a field line passing through 
y = 0, z = 0 may be readily obtained. For if a field line passes through y = 0, 
z = 0, at y we have z(y) = s(y) - s(0) = A (y). If we define the correlation function 
P (y) as 

< s (x ) s(x + y) > = < s(0) s(y)> s <s 2 > p (y) 


the probability density for s Q = s(0) and Sj = s(y) is given by 


f(s 0 , s^ds^Sj 


exp ^ (s * + sj - 2 s Q s j »/')/ <s 2 > ( 1 -</< 2 )| 


2n <s 2> [ 1 - i/< 2 ] 1 2 


ds Q ds j 


(13) 


where 


<sj> = <sj> = <s 2 > and p = p(.y) 


If we express s t as s t = s 0 + A in equation (13) we may then integrate over all 
values of s Q holding A fixed to obtain the probability density for A 



f (A)dA = 


exp - —• A 2 /< s 2 > (1 - «/<) 

L « ■ 

( 2 n <s 2 > (1 -»//)] 1/2 


dA 


(14) 


We see immediately that this is not the same as the probability density for e, 
under the condition that s 0 is zero. This would be given by 


f ( S J S 0> V 0 - [f(s | . s o)/f(s 0 )l , 0 s o 

exp [" T s i /<s2> (i -'A 2 )] 

[ 2 ^<s 2 > (1 - y 2 )] 1/2 

and can be seen to differ from expression (14) in containing the factor (1 -i' 2 ) 
rather than (1 - 0). The reason that expression (14) must be used rather than 
(15) is that the knowledge that a field line passes through the point y = 0, z = 0 
does not tell us that s = 0; we do not know the value of s Q and hence must 
average over all possible values. 

Our next problem is to find the probability that a field line has not reached a 
point z t above or below the central plane in a distance y along the mean field. 
Unfortunately the absorbing wall analogy will not serve us in this instance. 
There is a dynamic equation governing the development of the density in equa- 
tion (14). It is easy to verity that (14) is the Green function for the equation 


( A »y) 

By 


+ <s 2> ^' 


B 2 f (A, y ) 
BA 2 


(16) 
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with the initial condition f( \ 0) - b (A) and where = d«/'/dy. In the appendix 
we show that the diffusion equation (16) may be derived directly using the formal- 
ism of Kubo (1961). It can be seen that equation (16) is equivalent to a standard 
diffusion equation wit. diffusion coefficient <S2> and a "time" variable (1 - >p) 
which varies only between 0 and 1. This means that the diffusion only turns on 
for a finite time (actually length y in our case) of order L after which nothing 
more happens. In the absorbing wall problem a field line has a finite probability 
of reaching the wall in a distance of order L but if it doesn't reach the wall by 
this distance there is a good chance that it never will. 

The reason for this difficulty is that equation (16) describes a process that is 
decidedly non- Markovian; the evolution of the probability density f(A, y) depends 
not just on the immediate form of f but on how far f has developed from its initial 
form f(A, 0), e.g. on its past history. It can be seen that if one starts with an 
initial f(A, 0) =b (A) the solution relaxes asymptotically to the form of equation 
(14). However* this is not a stable solution for if one starts with this form for 
f(A, 0) equation (16) will relax it further to a Gaussian with a dispersion twice as 
large as before. So we see that the analogy with particle diffusion completely 
breaks down here and we must turn to another method for solving this problem. 

The answer is found in the solution of the zero crossing problem as given by Rice 

(1945). Rice shows that for a random function z(y) with a probability density for 
< t 

z and z = dz/dy at a given y, f(z, z ; y) the probability that the curve crosses 
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the zero axis with a positive slope between y and y + dy is given by 

dy I z' f (0. z' ; y) dz (17) 
*0 

Expression (17) may be immediately generalized to our case to determine the 
probability density for a field line displacement to cross the z = z, axis with 
positive slope (since the first crossing must be outward we are only interested 
in outward crossings). We have 


f (y ! * l ) dy = dy 



f(Zj . 


A' 


y) dA' . 


(18) 


The probability density f(z t ,A ' ; y) may be derived as follows: Starting with the 
trivariate Gaussian density for s 0 , e , and sj, 

f (s 0 . s , . s | : y) -r (2rt y 3 2 (M)' 1 ' 2 exp-^xRlT 1 x 

where 


- (Sg , S | | S | ) 


and 


< s 2 > 


1 1 <p ■/' 


l o 




(0) 


wo may express Sj = s 0 + A and s t - A'. We may then integrate over all possible 
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values of the 8 () in a straightforward though tedious calculation to obtain 


f(A,A';y) - f2« )“ 1 ( * S 2 -’ D )~ l exp 


* 2 ( 1 - y ) 

"7s? , D 


* j fe 

.V/*// - ' V|Q - 

< S^>D O t .2 


2 s 2 D 


( 19 ) 


where 


D *-2v, 0) (1 ~r) 


Expression (19) may now be inserted into expression (18) and the indicated 
integration may be performed analytically to obtain 


D ‘ 2 

Uy. - > .—.I— -" — J 


4-^(1 - i') 


~cxp 


tf^ 1 


4” 1 2 (1 - ip) 3 2 


( 20 ) 


l - erf I ^ 

,2(t 


where ^ 

e ••• z | / 2 s 2 > 12 ^ 

To proceed beyond this point an explicit form must be chosen for the correlation 
function •/> = </- (y)» This function is not known for the motions of the gas in our 
galaxy so we shall have to choose one as an ad hoc procedure. This should not 
be a serious disadvantage since we shaL see that the main features of the func- 
tion in expression (20) do not depend on the detailed form ofv'(y) but can be 
understood on the basis of the scale of the correlation function, or correlation 
length L. 
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If the correlation length is defined by 


L - j p(y ) dy 

*o 

we may choose a Gaussian as a typical form for >/• (y) to obtain 


<p(y) - exp (-^y 2 4L 2 ) = exp (-x 2 ). 


With this choice we may now write the probability (per unit x) as an explicit 
function of x 


f(x; S) 


exp (-48 2 D(x)) Dl 2(x) + xS exp(-x 2 ) exp [-8 2 /l-exp(-x 2 )] 
2v[l -exp(-X 2 )] ^1/2 [1 „ exp ( . x 2)] 3 2 


1 - erf 


-2x8 exp (-x 2 ) 


exp (-x 2 )) 12 D 1 2 (x). 


( 21 ) 


Where we have multiplied expression (20) by a factor of 2 to allow for the fact 
that we are interested in the two surfaces z - ± z • In Figure 1 we have plotted 
the function f(x, 8 ) for three different values of 8 . The choice 8 = 0.8 was made 
for purposes of comparison with the results of Jokipii and Parker (1969 b) for 
there is the same unconditional probability finding a given field line outside of 
this distance from the central plane as there was for finding it beyond one scale 
height (~ 130 pc) in their model. The major features of these curves are readily 
understoodt for values of x » 1 (y » L) the flat portion 
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f (x ; r ) > ( 2rr ) ' 1 exp ( -b 2 ) 


simply reflects the fact that initial conditions have become forgotten and the 
probability has relaxed to its asymptotic state. 

The peak at x ~B can be seen to be due entirely to the ballistic propagation of 
initial conditions. At y = 0 the slope A' has a probability density function 

f(A')dA' - f. 2rr<s 2 > <p" (0)]~ l 2 exp (A' 2 /2<s 2 >0" (0) ) dA'. (22) 


If we assume that the field line continues with the same slope A' until it interests 
the surface z = z, it will interest aty = y 1 3 /A' . Substituting A' = /y, 

into expression (22) we obtain a probability density for y^ using 0 " (0) « ^/2L 2 


f(yi)dy t 



exp (- z 1 2 L 2 /77y 1 2 <s 2 > )dy t 
U 2 <s 2 > /L 2 ) ,/2 


= 77 ” 1 2 (S/Xj 2 ) exp (-S 2 /Xj 2 ) dxj. 


(23) 


It is easy to see that this expression fits the initial part of the curves in Figure 
1, rather well at least up to the maximum at x = s . Of course as B gets larger 
tiie peak gets smaller because the ballistic propagation of initial conditions 
washes out as the field line has to go a distance of the order of L or further 
before interesting the plane z ■ z r 
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The probability expressed by the density f(x; & ), expression (21), is the proba- 
bility that the field line will cross the surfaces z - ± z, between x and x + dx 
with no concern as to whether or not it is the first crossing since x = 0. Since 
the probability that the field line has not crossed the surfaces anywhere between 
0 and x is given by 



the probability that a field line crosses the surfaces between x and x + dx for 
the first time is given by 


dF 

dx 


f (x; 5 ) exp 



f(x' 



(24) 


We call this expression dF/dx since the integral F is the distribution function 
for the "end point" X of the field line, i.e., 


P (X £ x) = F(x) = 1 - exp 



(25) 


dF/dx and F(x) are plotted in Figures 2 and 3 respectively. In Figure 3 we have 

marked the values of x for which F(x) = 1/2 for the three values of 8 . For 8 = 

0.8, the case corresponding to that considered by Jokipii and Parker, we see 

that x is slightly less than 4 or in other words the field line has a 50% chance 
1/2 

of getting out after about 4 correlation lengths. This is somewhat less than the 
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estimate of 6 correlation lengths made by Jokipii and Parker <1969 b). The 
discrepancy originates in the fact that the previous authors asked a somewhat 
different question, namely how many observations, each one a correlation length 
apart, must be made for there to be a 50% cliance of catching the field line out- 
side the surfaces z = ± z x . It is clear that this question does not count those 
cases where the field line goes out and comes back in between successive obser- 
vations and hence overestimates the distance required for the field line to pass 
through the critical surfaces for the first time. All in all the discrepancy does 
not appear to be a large one and it is a bit surprising that the results are as 
close as they are considering the rough nature of the original estimate. 

IV. CONCLUSION 

We have seen that the model of a Gaussian process for the galactic gas motions 
allows us to make exact calculations about the probability density of the galactic 
magnetic field. Whereas the model does not exactly represent reality the fact 
that no mathematical approximations need to be made assures us that our results 
are at least physically realizable and presumably not too far from the truth. 

With this approach we have been able to calculate the distribution function for 
the "end points" of the field line passing through the earth, end points meaning 
with respect to cosmic ray confinement. The point x 12 of 50% escape proba- 
bility was seen to be roughly consistent with a previous estimate. More important 
a probability density function was obtained which should prove useful in investi- 
gating the stochastic aspects of galactic cosmic rays. 
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APPENDIX 


We start with an initial density of field lines f(z, y) (for one field line f(z, y) = 
&(z - z 0 ) where z Q is the position of the field line). If the material carrying the 
field lines is displaced in the z direction as a function of position y as s(y) the 
density function becomes a function of y as 

f(z,y) = f 0 (z - s (y) + s(0)) (A-l) 

where f 0 { . 1 ) is the density function at y - 0. This may be written as a Taylor 
expansion using the Taylor expansion operator as 


f(z,y) - exp 


- (s(y) ~ s(0)) 


f 0 < z > 


(A-2) 


where expanding the exponential in a power series of its argument will generate 
the conventional Taylor series. 

We may now take an ensemble average of (A-2) over the ensemble of possible 
displacements sty) assuming that the initial f Q (z) is statistically independent 
from sty) we have 

<f(z.y)> =<expf- (s(y) - s (°»f“|> <f 0 (z) > (A-3) 


It has been shown (Kubo 1961) that if sty) represents a Gaussian process the 
average of the exponential operator greatly simplifies and we have 


(A-4) 


<Vxp 


-(s(y> - s(0) ) 



exp 


\ • (s(y) - s(0) ) 2 • 



assuming' s(y) > = 0. Expanding 


(s(y) - s(0)) 2 > = <s 2 (y) + s 2 (0) - 2s(y) s(0)> = 2<s , >(l-^(y) ) (A-5) 


so 


<f (z, y )> - exp 


<s 2 > (1 - yp) 





(A-6) 


If we differentiate equation (A-6) with respect to y we can obtain the differential 
equation for which (A-6) is the solution. 

<f 0 > r -<s 2 ><// — <f>.(A-7) 
dz 2 

We can immediately see that this is just equation (16) whose solution is expres- 
sion (14). If we had taken ensemble averages over only that subset of the ensemble 
for which s(0) = Owe would have obtained a more complicated differential equa- 
tion, one that is satisfied by expression (15) as we would expect. 


B<f > 2 di i> d 2 

. - - <S * > exp 

dy dy dz 2 


<S 2 > (1 


d 2 

dz 2 


25 


REFERENCES 


Davis, L. and Berge, C. L., 1968, in Stars and Stellar Systems, ed. G. P. Kuiper 
and B. M. Middlehurst, Vol. 7: Nebulae and Interstellar Matter (Chicago: 
University of Chicago Press), chapt. XV. 

Davis, L. and Greenstein, J. L., 1951, Ap. J., 114, 206. 

Feller, W., 1966, An Introduction to Probability Theory and its Applications, 

Vol. II (New York: John Wiley and Sons, Inc.), chapts. I and III. 

Hiltner, W. A., 1956, Ap. J. Suppl., 2, 389. 

Hulst, H. C. van de, 1967, Ann. Rev. Astr. and Ap., 5, 167. 

Jokipii, J. R., and Parker, E. N., 1969a, Ap. J., 155, 777. 

1969b, Ibid, 799 . 

Jones, F. C., 1970a, Acta Physica Hung. 29, Suppl. 1, 23. 

1970 b, Phys. Rev. Letters, 25, 1534. 

Kubo, R., 1961, in Fluctuation, Relaxation and Resonance in Magnetic Systems, 
ed. D. Ter Haar, (New York: Plenum Press), chapt. HI. 

Kulsrud, R. and Pearce, W., 1969, Ap. J., 156, 445. 

Parker, E. N«, 1965, Ap. J., 142, 584. 

26 


! * X 

ji 

i.H 


4 



4 


* 



'' % 


% 






1969, Ibid 157, 1129. 

Rice, S. O., 1944, Bell System Tech. J., 23, 282 
1945, Ibid, 2^, 46. 

Both of the papers by Rice are reprinted in Selected Papers on Noise and Stoch 
astic Processes, ed. N. Wax. (New York: Dover Publications, Inc. 1954). 

Roberts, W. W., Jr., and Yuan, C., 1970, Ap. J., 161 , 887. 

Spitzer, L., Jr., 1968, in Stars and Stellar Systems, ed. G. P. Kuiper ami B. M. 
Middlehurst, Vol. 7: Nebulae and Interstellar Matter, (Chicago: University of 
Chicago Press), chapt. I. 


FIGURE CAPTIONS 


Figure 1. Outward crossing rate per unit x Y/2L of field line crossing 
surfaces z = ± z,. Curves for three values of b - Zj/2 < s 2 >* 2 are shown 
as functions of x. 

Figure 2. Outward first crossing rate of field line. Probability density per unit 
x =v“ y/2L for three values of & = z j /2 < s 2 > 12 . 

Figure 3. Probability that first crossing point is less than x = 1 ' " y/2L for three 
values of b - z J 2 < s 2 >i /2 . 
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